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Abstract 

Based on the quantum molecular dynamics (QMD) picture, we cal- 
culate the energy and mass dependence of fragment production. For the 
present study, we simulated the reactions of 20 Ne + 20 Ne, 40 Ar + 45 Sc, 
58 Ni + 58 Ni, 86 Kr + 93 Nb, 129 Xe + 124 Sn and 197 Au + 197 Au at cen- 
tral geometry. Our findings clearly show a linear mass dependence for the 
peak center-of-mass energy at which the maximal IMF production occurs. 
Such linear dependence for peak center-of-mass energy on the system size 
has also been observed in recent experimental studies. We also predict 
a similar behavior for the multiplicities of different kinds of fragments. 
Experiments are called for to verify this prediction. 



Keywords: heavy-ion collisions, multifragmentation, quantum molecular dynam- 
ics (QMD) model. 

1 Introduction 

The central heavy-ion (HI) collisions have been reported to result into a complete 
disassembly of nuclear matter at bombarding energies above 100 MeV/nucleon. 
This disassembly of hot and dense nuclear matter also commonly known as the 
onset of multifragmentation, is found to occur when nuclear density drops to 
less than half of its initial value pQ. In the low density phase, the onset of 
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multifragmentation is expected to occur due to Coulomb instabilities. It is well 
known that the mean multiplicity of intermediate mass fragments (Nimf) de- 
pends strongly on the bombarding energy as well as on the impact parameter of 
the reaction. Leray et al. [T] studied the fragment distribution for the reaction 
of O+AgBr near the point of threshold multifragmentation and reported the 
onset of multifragmentation around 150 MeV/nucleon. Earlier, Peilert et al. [2\ 
have shown that the true multifragmentation events were confined to central 
collisions only. In their study, 93 Nb + 93 Nb reaction was simulated at incident 
energies between 30 and 200 MeV/nucleon. In central collisions of 93 Nb+ 93 Nb, 
one observes the maximal multiplicity (Mc(A > 4)) around 100 MeV/nucelon. 
More recently, an extensive and exhaustive study by Puri et al. [3j reported the 
outcome of 40 Ca + 40 Ca reaction over incident energies between 20 and 1000 
MeV/nucleon and over the entire impact parameter range. This study indi- 
cated the generation of events from incomplete fusion-fission to multifragment 
emission and finally complete disassembly of the nuclear matter. They observed 
a peak in the fragment production around 60 MeV/nucleon in the central col- 
lisions. Interestingly enough, a rise and fall in the multiplicity with impact 
parameter was not observed for low incident energies (20-40 MeV/nucleon). 
The existence of peak energy for maximal IMF emission was in accordance with 
an earlier study by Peilert et al. }2 using Au nuclei. The multifragmentation, 
therefore, exhibits a complex picture which is quite sensitive to the entrance 
channel characteristics i.e., to the impact parameter, beam energy as well as to 
the total mass of the target and projectile [TJ HI [5J [El [B] ■ The beam energy 
dependence of IMF emission was recently analyzed by Sisan et al. [9j using MSU 
47r-Array set up. In their study, emission of intermediate mass fragments was 
reported for the central collisions of 40 Ar + 45 Sc, 58 Ni + 58 Ni and 86 Kr + 93 Nb. 
They predicted a rise and fall in the emission of IMFs with beam energy and 
observed a linear peak energy dependence on the size of the system. The per- 
colation calculations used in the above study, however, could not fully explain 
this dependence. This led to the conclusion by Sisan et al. that perhaps phase 
space models can explain this observation [9] . We plan to address this situation 
by employing a dynamical model, where one can follow the reaction dynamics 
from the start to the end where matter is cold and fragmented. Our present 
study employs microscopic quantum molecular dynamics (QMD) model [101 111] 
which is described in section O Section [3] is devoted to model calculations and 
results, which are finally concluded in section 0J 
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2 Description of the model 



The quantum molecular dynamics (QMD) model is a time dependent A-body 
theory which is able to describe the many body phenomenon like fragment 
formation. Here each nucleon in phase space is represented by a Gaussian wave 
packet of the form: 



^i(r, Pi (i),ri(t)) = (2^3/4 ^ 
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(1) 



Mean position ri(t) and mean momentum Pi(t) are the two time dependent 
parameters. The Gaussian width \f~L is centered around the mean position 
rt(t) and mean momentum Pi(t) and is same for all nucleons. This value of \f~L 
corresponds to a root-mean-square radius of each nucleon. The effect of different 
Gaussian width in fragmentation is reported in reference |12| . The centroids 
of Gaussian wave packets (rf(t), Pi(t)) in phase space follow the Hamilton's 
equations of motion [TU1 HJ : 

or l dpi 

In the above equations, (H) stands for the total Hamiltonian of the system, 
which consists of kinetic and potential energy terms: 

A T +A P 2 -i A T +A P 

<*>= £ 1^ + 2 £ vir+v^ k +v^, (3) 

At and Ap being the target and projectile masses. The nucleon-nucleon in- 
teraction in consists of a local Skyrme interaction, a long-range Yukawa 
interaction and an effective charge Coulomb interaction parts [101 lllj : 

V}° c = t 1 5(r l -rj)+t 2 5(r t -r ] )5(r l -r k ) (4) 
yYuk = t «sp{-k<-*j|}/A* 
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Zi, Zj are the effective charge of baryons i and j. In QMD model, one neglects 
the isospin dependence of the interaction. All nucleons in a nucleus are assigned 
the effective charge Z = jf^rg^ jTT] • The long-range Yukawa force is necessary 
to improve the surface properties of the interaction. The parameters /i, t\,t2, £3 
in ((4]) are adjusted and fitted so as to achieve the correct binding energy and 
mean square root values of the radius of the nucleus [10]. Since QMD model 
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follows the time evolution of nucleons only, one has to construct the fragments. 
In a simplest approach, two nucleons are assumed to share the same cluster if 
they are closer than a distance of 4 fm. This method, also known as minimum 
spanning tree (MST), can be applied when matter is dilute and well separated. 
This picture is true when incident energy is high and collisions are central in 
nature. One has to also keep in the mind that semi-classical models like QMD 
can not keep nuclei stable for long time. A typical stability of nuclei can be 
seen untill 200 fm/c. If one analyzes the fragment formation with MST alone, 
then one may not achieve true fragment structure at 200 fm/c. To speed up the 
recognition of fragment structure, we add secondary condition that fragments 
produced with MST method are subjected to further binding energy check: 



i=l 



( Pi - P%?) 2 1^„ , s 



2mi 2 



< Eund- (5) 



We take Eund — -4.0 MeV/nucleon if > 3 and E^nd = otherwise. In 
this equation, N' is the number of nucleons in a fragment, PS^ is the center-of- 
mass momentum of the fragment. This modified version of conventional MST 
method with binding energy check is labeled as MSTB method. The magnitude 
-4.0 MeV/nucleon of Ebi n d is able to recognize the fragment structure quite ac- 
curately. It is chosen keeping in the mind the average binding energy of clusters. 
In a recent communication 13J, we used instead microscopic binding energies 
based on experimental information. Nearly no effect was seen by varying the 
binding energy. We have shown in many calculations that this check is close to 
other momentum cuts [TH [THl [T7] or sophisticated algorithms like simulated 
annealing clusterization algorithm [T51 Qjjl [20] . 

We employ a soft equation of state (EoS) along with Cugnon parametrization 
of n-n cross section for the present study fTUj. The choice of soft EoS has been 
advocated in many theoretical studies. Recently, Magestro et al. [21] tried to 
pin down the nuclear incompressibility using balance energy. Their detailed 
study pointed towards a softer equation of state. Another study concerning the 
linear momentum transfer occuring in central HI collisions [22] also showed that 
a soft compressibility modulus is needed to explain the experimental data. 



3 Results and Discussion 

Here, we simulate the central heavy-ion collisions of 20 Ae + 20 Ae (Ei a b=10 to 
55 AMeV), 40 Ar + 45 Sc (E lab =35 to 115 AMeV), 58 Ni + 58 Ni (E lab =35 to 95 
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AMeV) and m Kr + 93 Nb (Ei ab =35 to 95 AMeV), 129 Xe + 124 Sn {Ei ab =45 to 
130 AMeV) and 197 Au + 197 Au (E lab =70 to 130 AMeV). The systematic study 
over a wide range of beam energies and system masses allows one to confront 
the theoretical predictions with experimental findings and search for the mass 
dependence. Note that only symmetric reactions are taken for present analysis. 
Our calculations are performed at fixed impact parameter b= fm. We calculate 
the reaction at fixed incident energies and then calculate corresponding center 
of mass energy. For each such set, 500 events were simulated that minimizes the 
fluctuations to greater extent. The choice of central collisions for the present 
study guarantees the formation of highly excited systems that may break into 
a large number of pieces. Further, the emission from such events is almost 
isotropic, which may represent a 'single source' emission. In figure [TJ we display 
the average multiplicity of intermediate mass fragments (Nimf) calculated as 
a function of beam energy £ c , m . in the center-of-mass frame employing MSTB 
method. We display here the model calculations for unfiltered events of four 
entrance channels 40 Ar + 45 Sc, 58 Ni + 58 Ni, 86 Kr + 93 Nb and 129 Ae + 124 Sn. 
The {Njmf) first increases with beam energy, reaches a peak value and then 
decreases. This trend is visible in all of the four entrance channels shown here. 
This trend is less clear for the lighter 4:0 Ar + 45 Sc system whereas it is more 
clearly visible for the heavier systems. A similar dependence of (Nimf) on 
center of mass energy is also observed in experimental data taken with the MSU 
47r- Array [9] . This behavior can be understood in terms of compression energy 
of the system. With the rise in the beam energy, compression energy breaks 
the IMFs into lighter mass fragments thereby, leading to fall in the multiplicity 
of IMFs. The maximal E cm _ and corresponding peak (Nimf) was obtained 
through a quadratic fit to the model calculations. One should also note that 
the shape of the beam energy dependence of IMF production is quite close to 
one reported in the experimental data [9]. As reported by Sisan et al. [9], 
the peak E c . m . extracted for different entrance channels scales with the size of 
the system. Such scaling is also visible in our present calculations (see Fig. [1}. 
In figure [2) we display the time evolution of average density along with the 
mean multiplicity of intermediate mass fragments (Nimf) defined as fragments 
3 < Z < 20 plotted at the peak center-of-mass energy E c _ m , (at which maximal 
IMF emission occurs). We now include 20 Ne+ 20 Ne and 197 Au+ 197 Au systems 
also for the study of system size effects. The average nucleonic density of the 
system is calculated as: 
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Figure 1: The mean IMF multiplicity (Njmf) versus beam energy E c . m . for the 
reaction of 40 Ar + 45 S*c, 58 Ni + 58 Ni, 86 Kr + 93 Nb and 129 Xe + 124 Sn. Open 
circles depict the calculations employing QMD + MSTB approach for unfiltered 
events. The quadratic fits (solid curves) to the model calculations are drawn to 
estimate the peak energy at which the maximal IMF emission occurs. 
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Figure 2: The time evolution of mean nucleon density (upper panel) and mean 
IMF multiplicity (lower panel). Results displayed here are at the energy for 
peak IMF production. 
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Figure 3: The system size dependence of the peak E c . m . and peak (Nimf)- 
Our model calculations (open circles) for unfiltered events are compared with 
experimental data (solid squares). Also shown in the figure are the percolation 
calculations (open squares) [9]. 



with Ti and Tj being the position coordinates of the i th and j th nucleons. 
The Gaussian width L is fixed with a standard value of 1.08 fm 2 . As expected, 
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the average nucleonic density has a mass dependence, being maximal for the 
197 Au + 197 Au system and minimal for the 20 Ne + 20 Ne system. This also 
indicates a linear density dependence on the system size. The intermediate 
mass fragments also show similar mass dependence. One can also notice that 
fragment production almost saturates around 200 fm/c. In other words, time 
span of 200 fm/c is large enough to pin down the fragment structure. The 
maximal fragment production is for 197 ^4u+ 197 ^4u system whereas 20 Ne+ 20 Ne 
system results in minimum value. It may be mentioned that IMF multiplicities 
obtained in 20 Ne+ 20 Ne and 40 Ar+ 45 Sc collisions exclude the largest and second 
largest fragment respectively to infer the system size dependence accurately. 

We plot in figure 02 the peak E c , m , as well as peak (Njmf) as a function 
of total mass of the system Atot- Strikingly, our model calculations employ- 
ing MSTB approach are in good agreement with the experimental data (solid 
squares) of MSU 47r-Array for peak E c , m _. For peak (Nimf), some deviation 
can be seen for heavier masses. This could also be due to the fact that our 
calculations are not filtered for experimental acceptance. 

One can also see that the predictions of percolation model fail to explain 
the sharp dependence of peak E c _ m . on system mass. Our present results show 
a linear mass dependence of the form: mA to t + c for the peak E c , m _. These 
observations suggest that the peak E c _ m ., thus, acts as a measure of finite size 
effect. It is worth mentioning that the critical excitation energy was estimated 
from the cluster size distribution fitted to power law: u(A) cx A~ x at different 
beam energies for which the exponent A reaches a minimum. Based on the 
percolation calculations, the critical excitation energy is also found to increase 
when initial lattice size increases [TJ. Interestingly, the mass scaling of peak 
(Njmf) can be reproduced with a power law: cA[ ot with exponent close to 
unity. In figure 31 we finally extend the above study for various fragments 
consisting of free nucleons, fragments with mass A=2, light charged particles 
LCPs [2 < A < 4] , medium mass fragments MMFs [5 < A < 9] as well as heavy 
mass fragments HMFs [10 < A < 44]. Interestingly, in all the above cases, 
a clear system size dependence can be seen in a manner similar to (Nimf) 
dependence. We observe a power law of the form cA T tot ; A to t is mass of the 
composite. In all the cases, parameter r is very close to unity. As noted in 
[5] , the percolation model failed badly to reproduce the power law dependence. 
A linear mass dependence observed with value of r ~ 1 depicts the picture of 
vanishing surface- Coulomb effects. Experiments are called for to verify this new 
prediction. 
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Figure 4: The multiplicity of (a) free nucleons, (b) fragments with mass A=2, 
(c) light charge particles LCPs, (d) medium mass fragments MMFs, and (e) 
heavy mass fragments HMFs as a function of total mass of the system A tot . 
Model calculations done at peak E c , m , (open circles) are fitted with power law 
of the form: cAJ ot . 



4 Summary 

We aimed to reveal the dependence of IMF production on beam energy and 
system size. This was achieved by a study over wide range of system masses 
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and incident energies. Our present results reproduce the experimental trend 
of both rise and fall in (Njmf) with beam energy. At the point of onset of 
multifragmentation, we obtained the scaling of peak E c , m . with system mass. 
The observed trend of peak center of mass energy is in agreement with previous 
experimental studies [7J [221 (Ml (HI The behavior of average nucleon 

density as well as IMF multiplicity at the point of onset of multifragmenta- 
tion for different entrance channels also reflect the dominance of system size 
effects. Our calculations also reproduce a power law of the form cA[ ot ; A to t 
being the total mass of the system. We predict a similar power law dependence 
for the fragments of different sizes at the energy for peak IMF emission. Inter- 
estingly, as observed experimentally, the exponent r is close to unity in all cases. 

This work is supported by CSIR, Government of India vide grant no. 7167/NS- 
EMR-II/2006, India. 
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